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PG 1613+426: a new sdB pulsator* 
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Abstract. We report the detection of short period oscillations in the hot subdwarf B (sdB) star PG 1613+426 from 
time-series photometry carried out with the 91-cm Cassegrain telescope of the Catania Astrophysical Observatory. 
This star, which is brighter than the average of the presently known sdB pulsators, with B = 14.14 mag, has 
Tcff = 34 400 K and logy = 5.97, its position is near the hot end of the sdB instability strip, and it is a pulsator 
with a well observed peak in the power spectrum at 144.18 ± 0.06 s. This star seems to be well suited for high 
precision measurements, which could detect a possible multi-mode pulsation behaviour. 
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1. Introduction 

The hot subdwarf B (sdB) stars form a homogeneous 
group populating an extension of the extreme horizontal 
branch (EHB) in the (T c g — log (^-diagram towards tem- 
peratures up to 40 000 K. These stars are evolved low-mass 
(~ 0.5M Q ) objects with a He burning core surrounded by 
a H surface layer which is too thin (< 4% by mass) to sus- 
tain the H-shcll burning. Their origin is still debated, but 
it seems that they have experienced He-flash phase and 
a substantial mass loss along the red giant branch. Their 
further evolution proceeds toward the EHB by crossing 
the subdwarf O population and eventually entering the 
white-dwarf graveyard (Maxted et al. 2001, Heber et al. 
2002). 

The recent discovery that several of them are multi- 
mode pulsators has triggered both observational and the- 
oretical efforts for studying their characteristics and pul- 
sating mechanisms (Charpinet et. al. 2001). Data on 30 
of such pulsators, also known as sdB variables (sdBVs), 
are reported by Charpinet (2001), including the sdBV 
HS0702+6043, the detailed results of which are reported 
in Dreizler et al. (2002), and the two sdBVs PG1325+101 
and PG2303+019, recently discovered by Silvotti et al. 
(2002). The 31st star of this class was discovered by 
Piccioni et al. (2000) and later confirmed by Ulla et al. 
(2001). These stars show short period (1 — 10 min) and low 
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amplitude (1-50 milli-magnitude [mma]) non-radial pulsa- 
tion modes. Charpinet et al. (1996,1997) have shown that 
the non-radial pulsations are probably driven by an opac- 
ity bump associated with iron ionization. 

Here we report the discovery of a new variable of 
this class, the 32nd one, namely the PG 1613+426 sdB, 
thanks to the data collected at the Catania Astrophysical 
Observatory. This star, with B = 14.14 mag, is brighter 
than the average of the presently known sdB pulsators. 
This object has been selected from the list of ultraviolet- 
excess Palomar Green objects (Saffer et. al 1994), where 
a determination of the main atmospheric parameters 
is available. For this star, low resolution spectroscopy 
(~ 6 A) leads to a sdOA spectral classification due to its 
strong, broad Balmer and Hel absorption lines, and in par- 
ticular gives T off = 34400 ± 500 K, Xogg = 5.97+0.12, and 
N(He)/N(H) = 0.022 ± 0.003 from a least-square fitting 
of Balmer lines with LTE models. These characteristics 
place this objects near the hot end of the theoretical sd- 
BVs instability strip. 

2. Photometric observations 

The photometric observations were carried out with the 
the 91-cm Cassegrain telescope of the M. G. Fracastoro 
stellar station of the Catania Astrophysical Observatory 
(Serra la Nave, Mt. Etna, 1720 m a.s.l.). The telescope was 
equipped with a photon counting photometer, which used 
an unfiltered EMI 9789-QA photomultiplier as detector. 
A 21" diaphragm was used to isolate the star light from 
the sky background. Owing to the PG 1613+426 spectral 
distribution and the response curve of the photocathode, 
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Table 1. Summary of the observations. The initials of the 
observers reported below are: AB = A. Bonanno, AF = 
A. Frasca, GM = G. Mignemi. 
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the photometric passband included from B to UV wave- 
lengths. Data were collected by integrating in time inter- 
vals of 15 or 20 s, obtaining typical counts of about 20 000 
per measurement. The estimated photometric accuracy at 
our observing station was approximately 5 mma in good 
observing nights. Exposures of a reference star in the same 
stellar field and sky background were also taken during the 
night in order to subtract the sky background and check 
the overall stability of the photometric system. The cor- 
rection for atmospheric extinction was applied by using 
the average seasonal extinction coefficient for B Johnson 
filter. Since our spectral band did not coincide with the 
Johnson B band and the variation of the sky transparency 
was not fully taken into account by the non-simultaneous 
measurements of the comparison star, possible long-term 
trends of differential magnitude were removed by subtract- 
ing a low-order polynomial curve fitted to the data of the 
individual runs. 



3. Data analysis and results 

PG 1613+426 was identified as a variable star on 2002 
June 14. We have thus decided to observe again this tar- 
get on June 15 and June 16. Further observing runs were 
then performed in August 8, 9, 13, 14 and 15 in order to 
resolve the amplitude spectra with a better accuracy. A 
summary of the observations, including the starting time, 
date, duration, exposure time per measurement, name of 
the observers and the main pulsation period, is reported 
in Tabic 1, while a sample of the differential light curves 
is shown in Fig. [l]. 

The results concerning the amplitude spectra are il- 
lustrated in Fig. |^ for the three individual August 13, 
14 and 15 nights and the June 16 night. The amplitude 
spectrum of the three August consecutive nights, during 
which the observing conditions were particularly good, 
taken toghether, is shown in Fig. |^. The amplitude spec- 
trum shows a clear, dominant peak at 6.936 ± 0.003 mHz 
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Fig. 1. Sample of the differential light curves of 
PG 1613+426. Long term trends have been subtracted 
with low order polynomial fitting. 



(144.18 ± 0.06 s) with an amplitude of about 5 mma in 
each observation night. The period corresponding to the 
highest peak in the power spectrum is reported in Table 
1 for each night. The error has been estimated from the 
FWHM of the main peak in the spectral wind ow. The 
Lomb-Scargle periodogram analysis (Scargle 1982 ) applied 
to our data gives a confidence level > 99.9999 % for the 
highest peak in the spectrum obtained in the June 16, and 
August 13, 14 and 15 observation runs. 

We have refined our analysis by pre-withening the light 
curves with a sinusoidal fit of the dominant period, and 
iterating to obtain the light curve residuals for any new 
frequency which was found []. We found additional fre- 
quencies, but we noticed that the associated reduced \ 2 
did not change significantly during the pre-whitening pro- 
cedure, indicating that, at the present level of signal to 
noise ratio, the analysis precludes any further peak iden- 
tification. 

As a complementary method we used the CLEAN it- 
erative deconvolution algorithm (Roberts et al. 



1986) in 



order to eliminate the effects of the observational spectral 



1 We used the software Period98 by M. Sperl, University of 
Vienna, 1998. 
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Fig. 2. Amplitude spectra for some observing nights. The 
spectral window is shown in the insets. 
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Fig. 3. Amplitude spectrum for the three August consec- 
utive nights. The spectral window and a detail of the spec- 
trum around its maximum are shown in the insets. 

window in the power spectrum. The cleaned periodograms 
of the individual nights do not differ significantly from the 
non-cleaned spectra, since the spectral windows for essen- 
tially equally spaced data are always rather clean with 
side-lobes of low amplitude. The results of this analysis, 
applied to the three August consecutive nights, do not re- 
veal any difference from the previous one, confirming the 
clear presence of the peak at 144.18 mHz. 

Although the accuracy of our measurements is too 
low to reach a firm conclusion on the existence of ad- 
ditional periods, we cannot rule out the possibility that 



PG 1613+426 is a multimode pulsator. For instance, we 
note that frequencies in the range 6.63 — 7.92 mHz are ex- 
pected on the basis of the theoretical models discussed in 
Charpinet et al. (2001) and on observations of very simi- 
lar spectra in other hot sdBVs which have nearly the same 
position in the sdBVs instability strip, like HS 0815+4243, 
with T cS = 33 700 K and log. 9 = 5.95, and HS 2149+0847, 
with T off = 35 600 K and logs = 5.9 (0stensen et al. 
2001). 

We think that it would be important to observe this 
object by means of high precision CCD photometry in 
order to resolve its spectrum with better accuracy, and 
we plan to pursue this project with observations at our 
site in the near future. 
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